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1. Introduction 

The understanding of the vacuum structure of QCD is an important part of theoretical particle 
physics. The distribution of the topological charge density is believed to describe nonperturbative 
phenomena like the large r\' mass through the axial U(l) anomaly and the spontaneous breaking 
of chiral symmetry. The relation to confinement is less obvious. The instanton picture models the 
QCD vacuum as a random (quenched) or interacting (full QCD) instanton-antiinstanton system, 
whereas the understanding of confinement seems to require singular fields emanating from lower- 
dimensional defects [pp. 

Recent advances in implementing chiral symmetry on the lattice have made it possible to find, 
from first principles, manifestations of these dual aspects of vacuum structure. Overlap fermions 
are described by the Neuberger-Dirac operator constructed as a solution of the Ginsparg-Wilson 
relation. They are the cleanest known implementation of lattice fermions. For gauge fields free of 
dislocations, the index theorem is unambiguously realized and the spectrum consists of chiral zero 
modes and non-chiral non-zero modes. While the number and chirality of zero modes provides a 
value for the total topological charge, the non-chiral non-zero modes appear in pairs of eigenvalues 
±iX. 

For our methodical study we use quenched configurations simulated by means of the Liischer- 
Weisz action. In the following table we list the statistics of lattices used in our investigation. We 
show also the topological susceptibility % t . To set the scale, we use the quenched pion decay 
constant/^ = 91(2) MeV. 



p 


a [fm] 
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250 
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0.138(3) 


12 3 x 24 


1.66 


15.04 


0.647(58) 


251 


137-154 


8.00 


0.154(4) 


16 3 x 32 


2.46 


73.72 


0.607(21) 


2156 


160-178 



After presenting results concerning localization and chiral properties of individual eigenmodes 
in different parts of the spectrum up to A = 400 MeV, we describe two specific applications of the 
overlap Dirac operator. One is aiming to extract a UV-filtered field strength tensor [Q] in order to 
assess the local degree of (anti-)selfduality. The other uses the fermionic definition of topological 
charge density [g]. We interpret the respective structures and their dimensionality. 

2. Localization, dimensionality and local chirality of the lowest eigenmodes 

Recently, the localization, the eventually lower dimensionality and the local chirality of the 
lowest eigenmodes of the Dirac operator have attracted a lot of interest [|], ||]. The reason for this 
activity is that these modes might be pinned down by some topological defects. This is inspired by 
certain confinement mechanisms, e.g. the center vortex mechanism. In model field configurations 
of thick vortices the zero modes have been seen [|6|] localized along the vortices with additional 
peaks at the intersections. The particular interest in the lowest eigenmodes is motivated by the 



2 



The QCD vacuum probed by overlap fe unions 



Volker Weinberg 



observation that the lightest hadron propagators are very well approximated by taking only 0(50) 
lowest eigenmodes in the quark propagator into account. 

A useful measure to quantify the localization of eigenmodes is the inverse participation ratio 
(IPR), / = V Y, x P ( x ) 2 > w i tn the scalar density p (x) = \\r\ (x) \jfx (x) for eigenf unctions with eigenval- 
ues iX, which are conventionally normalized as Y*xP{ x ) = 1- In Fig. [I] (a) we plot the IPR averaged 
over bins with a bin width AA = 50 MeV and for the zero modes considered separately. The av- 
erage IPR shows an L and a dependence for zero modes and for non-zero modes only in the range 
A < 150 MeV. Beyond that the average IPR corresponds to 4-dimensionally extended states. In the 
theory of the metal-insulator transition, a more detailed description of the localization properties at 
different "heights" of the wave functions has led to consider generalized IPRs like I p = Y,xP( x ) p > 
averaged over parts of the spectrum and different realizations of disorder. In any physical appli- 
cation, dimensionalities d* (p) different from the embedding dimension d and eventually varying 
with p, can be inferred from the volume scaling, 



I p oc l-^p-^ 



L -d'{ P ){p-\) 



I p oc const 



extended wave functions, 
multifractal wave functions, 
localized wave functions. 



metallic phase, 
critical region, 
insulator phase, 

While the dimension for the metallic or insulator phase is d or 0, respectively, the critical multi- 
fractal region in between is characterized by a multitude of fractal dimensions d*(p) and also by 
critical level statistics ]M. 
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Figure 1: (a) The average IPR is shown for zero modes and for non-zero modes in X bins of width 50 MeV 
for the five ensembles, (b) The fractal dimensionality d* (p) obtained from fits of the volume dependence of 
the averages of the generalized /„ is presented for zero modes and for non-zero modes in X bins of width 
50 MeV for the three ensembles with different volumes and common j3 =8.45. 

From a fit of the volume dependence of I p we find that all eigenstates with X < 200 MeV 
in this sense belong to a critical regime. From the standard IPR, I = V h, we conclude that the 
dimensionality of the zero modes is close to two. In the next bins of width 50 MeV each, the 
dimensionality is near three, in agreement with [Q] and disagreement with [j|] . The IPR of the 
higher modes is independent of V and a, describing the observation that the modes freely extend 
throughout d*(2) < d = 4 dimensions. 

'However, in [Q] the Asqtad Dirac operator on quenched Symanzik improved gauge fields is used, while the analysis 
in [fin is done with overlap fermions on quenched SU(2) configurations generated with the standard Wilson action. 
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Looking at the higher moments of IPR one can see a decrease of the multifractal dimension 
with increasing p as shown in Fig. [I] (b). This reflects the fact that I p with p > 2 explores regions 
of higher scalar density. For p = 20 . . . 30 the emerging fractal dimension seems to converge. The 
envelope reveals that the maxima of the scalar density of zero modes and the lowest non-zero 
modes up to X = 100 MeV are characterized by d*(p) < 1 localization (pointlike peaks), whereas 
in the spectral region 100 MeV < X < 400 MeV, d*(p) = 2 or 3 characterizes the maxima of the 
scalar density. In contrast to this, the standard IPR alone would describe the modes in the spectral 
region X > 150 MeV as four-dimensional. For the zero modes, the results favor an interpretation 
in the vortex picture above, but that needs direct confirmation. 




(a) (b) (c) (d) 

Figure 2: Normalized histograms of the local chirality X(x) in the Q = subsample (consisting of 37 
configurations) among the 16 3 x 32 lattices generated at j3 = 8.45. (a) and (c): for the lowest ten pairs 
(with positive and negative X) averaged over the subsample. (b) and (d): for all modes from the subsample, 
averaged over eight X bins of width 100 MeV. (a) and (b): with a cut for 1/100 of the sites with biggest 
scalar density p (x). (c) and (d): with a cut for 1 /2 of the sites. 

The local chiral density pxs(x) = V / ^( x )7sV / a( x ) with L.tPasC*) = ±1 or for chiral or non- 
chiral modes, respectively, is equally important. Whereas the zero modes are entirely chiral, the 
non-zero modes with globally vanishing chirality may still have a rich local chirality structure cor- 
related with the local (anti-)selfduality of the background gauge field. A measure proposed by Hor- 
vath [||] to describe the local chirality Pi±(x) = \j/l(x)P±Yx( x ) (with projectors P± = (1 ± 75) /2 
onto positive/negative chirality) is given by 

X(x) = — arctan (r% (x)) — 1 £ [-l,+ll , (2.1) 

71 

with rxi x ) = Pk+{ x ) I 'Pa-M- F° r the chiral zero modes X(x) = ±1, while for the non-chiral non- 
zero modes X(x) is a strongly varying field. 

In Fig. H normalized histograms with respect to X (x) are shown for the Q = subsample (37 
configurations) of the 16 3 x 32 lattices generated at j3 = 8.45, with cuts including 1/100 and 1/2 
of the sites according to the biggest p (x). For zero modes, the histograms would be concentrated at 
X = ±1 irrespective of any cut. For the non-zero modes, it turns out that with a more restrictive cut 
in p(jc) the local chirality, i.e. the concentration of the histogram over the included lattice points, 
can be driven towards X = ±1. A similar enhancement can be seen towards the lower modes. For 
the lowest two pairs of non-zero modes all sites with p(x) above median show the two chiral peaks. 
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3. The UV-filtered gluonic field strength tensor and the topological charge density 

The lowest eigenmodes of the overlap operator can directly be used as a filter of the underlying 
gauge fields. Following a method proposed by Gattringer [0], the gluonic field strength tensor can 
be represented as 

^vW^I^/^W). (3-1) 

j 

The contribution of the j-th eigenmode f^ v (x)j = — | ^{x)y^y v T a \ifj{x) to the field strength tensor 
F^ v {x) is projected out by means of an appropriate combination of y-matrices with the SU(3) 
generator T a . To study the local (anti-)selfduality of the filtered field strength tensor, we analyze 
the ratio 

r(x) = {s(x)-q(x))/(s{x)+q(x)) , (3.2) 

with the action density ${x) = TrF llv {x) F^ v (x) and the topological density q(x) = TrF^ v (jc) F^ v {x) 
of the filtered field strength summing in Eq. ( |3.1| ) over different sets of low-lying non-zero modes. 
Using the ratio r(x) we do not need to know the normalization factor of s(x) and q(x). 



Mapping r(x) toR(x) E [—1,-1-1] analogously toEq. (2T), one would get R(x) = ±1 for totally 
(anti-)selfdual fields. In Fig. |we present histograms with respect to R(x) applying two different 
cuts with respect to the filtered action density s(x). We observe peaks at R = ±1 that become 
weaker with the inclusion of more modes. The peaks become more pronounced again when one 
concentrates on lattice points above some cut with respect to the filtered action density s(x). 




Figure 3: Normalized histograms with respect to the local (anti-)selfduality R(x) of the field strength tensor 
in the Q — subsample of 37 configurations generated on the 16 3 x 32 lattice at j3 = 8.45, (a) taken over 
the full volume, but depending on the number of non-zero modes (2 - 20) included in the filter, (b) applying 
various cuts according to the action density for a fixed number (20) of eigenmodes included in the filter. 

Alternatively to q{x) based on the filtered field strength tensor we also examine the topologi- 
cal charge density as given by the formula q(x) = — Tr [75 (l — | D;v(x,x))] [^] which holds for 
any 75-Hermitean Dirac operator satisfying the Ginsparg-Wilson relation. We use it in a two- 
fold way [TO) ], (i) We directly calculate the trace of the overlap operator. The "all-scale" 
density q{x) computed in this numerically very expensive way includes charge fluctuations at 
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all scales down to the lattice spacing a. (ii) We use the mode-truncated spectral representation 

«u«(*) = -Ew<;u(i ~ h) PasW- 

In Ref. [11] we have begun a cluster analysis of the two forms of the topological density. The 
number of separable clusters of q(x) grows without limit with a — > 0, and similarly the number 
of clusters of qi with increasing A cu t at fixed a. We have also seen [ |l2| ] that the negativity of 
the two-point-function of the topological charge density fll3| ] (up to a positive core due to the non- 
ultralocality of q{x)) is realized for the all-scale density only and increases with decreasing lattice 
spacing a. The correlator for q* describes clusters of the mode-truncated density. For sufficiently 
large Ac Ut , clusters of different sign are attractively correlated as demonstrated by a negative tail of 
the corresponding correlation function. 

In Fig. |] (a)-(c) we compare in a two dimensional cut the density q^ M with Ac Ut = 200 MeV, 
the density q(x) from the filtered field strength tensor and the all-scale density q(x). Viewed is a 
typical configuration generated at j3 = 8.45 on a 16 x 32 lattice. Normalizing q{x) in such a way 
that the maxima of q(x) and q^ix) are equal, we see that the two definitions (both including 18 
modes) are in very good agreement. The density of clusters is 0(1 fm~ 4 ), and the number of modes 
necessary to see the clusters approximately corresponds to the number of clusters. The topological 
susceptibility is smaller by a factor of 1/2 compared to the dilute instanton-gas approximation. 
This is an independent indication for pairing correlations among oppositely charged clusters of 

In contrast the cut through the all-scale topological density seems to be completely random. 
It is interesting to notice that structure can be extracted from this "noise" as well [14, 11]. First 
of all, cluster percolation at low q(x) leads to the two thin (thickness of order 0(a)) sign-coherent 
global clusters for the first time seen by Horvath et al. [JTOQ. At higher |g(x)|, exploring the interior 
of the clusters leads to the result that the topological landscape is also multifractal, with fractal 
dimensions ranging from d* = 0.7 to d* = 2.6. 

Finally, Fig. |] (d) and (e) show the positive and negative part of the local (anti-)selfduality 
variable R(x). One can see that the peaks of the mode-truncated topological density coincide with 
the domains of (anti-)selfduality. This is evidence that with the chosen smearing parameters, ex- 
tended lumps of high topological charge density may be interpreted as (anti-)selfdual objects. On 
the other hand, the all-scale density allows one simultaneously to localize singular fields with frac- 
tal dimension d* < 4. A more quantitative analysis based on our cluster algorithm as drafted in 
flTlj] will be discussed in detail in a forthcoming publication. 
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(a) 



(b) 



(c) 




(d) 



(e) 



Figure 4: Two-dimensional cuts through a typical configuration generated at j3 = 8.45 on a 16 3 x 32 lattice. 
Direct comparison of (a) q{x) based on the filtered field strength tensor, (b) qx^Jx) as the mode-truncated 
density and (c) the all-scale charge density q(x). Subpanels (d) and (e) show the positive and negative part 
of the local (anti-)selfduality variable R(x). Except for (c) always 18 eigenmodes are taken into account. 
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